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Abstract 

Many abelian gauge theories in three dimensions flow to interact- 
ing conformal field theories in the infrared. We define a new class of 
local operators in these conformal field theories which are not polyno- 
mial in the fundamental fields and create topological disorder. They 
can be regarded as higher-dimensional analogues of twist and winding- 
state operators in free 2d CFTs. We call them monopole operators 
for reasons explained in the text. The importance of monopole opera- 
tors is that in the Higgs phase, they create Abrikosov-Nielsen-Olesen 
vortices. We study properties of these operators in three-dimensional 
QED using large Nf expansion. In particular, we show that monopole 
operators belong to representations of the conformal group whose pri- 
maries have dimension of order Nf. We also show that monopole 
operators transform non-trivially under the flavor symmetry group, 
with the precise representation depending on the value of the Chern- 
Simons coupling. 
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1 Introduction 



One of the most fascinating problems in quantum field theory is understand- 
ing non-perturbative equivalences ( "dualities" ) between superficially very dif- 
ferent theories. A classic example is the quantum equivalence of the massive 
Thirring and sine-Gordon models |2j. More recently, a number of dualities 
has been conjectured for supersymmetric gauge theories in three and four 
dimensions. The earliest proposal of this sort is the S-duality of iV = 4 d = 4 
super- Yang- Mills theory || £|, ||] . A decade and a half later, N. Seiberg pro- 
posed a dual description for the 4d CFT which arises as the infrared limit of 
N — 1 d — 4 super-QCD M. The dual theory is again the infrared limit of an 
N — 1 d — 4 gauge theory. This proposal generated tremendous excitement, 
and soon many other candidate dualities have been found (see Refs. |7|, |3| 
for a review). Later it was realized that many field-theoretic dualities follow 
from string theory dualities. 

Until now, all dualities in dimensions higher than two remain conjectural, 
and the physical reason for their existence is not completely understood. On 
the other hand, 2d dualities have a rather transparent physical meaning. 
For example, the sine-Gordon model has topological solitons (kinks), and it 
can be shown that a certain local operator which creates a kink satisfies the 
equations of motion of the massive Thirring model 0. It is believed that 
many higher- dimensional dualities arise in a similar manner, by "rewriting" 
the theory in terms of operators which create topological disorder. But it 
proved very hard to make this idea precise. 

There are several related difficulties that one encounters in dimension 
higher than two. First of all, interesting higher- dimensional dualities involve 
gauge theories. This implies that in order to write down an operator describ- 
ing the dual degrees of freedom, one has to work in an enlarged state space 
which includes the unphysical degrees of freedom of both the original and the 
dual gauge fields. It is not known how to construct such an enlarged space. 
Fortunately, there are non-trivial examples of dualities in three dimensions 
(so called 3d mirror pairs |9]]) for some of which the dual theory has a trivial 
gauge group. In this case one can hope to construct the operators describing 
the dual degrees of freedom directly in the state space of the original gauge 
theory. 

The second difficulty is that it is hard to construct topological disorder 
operators in interacting fields theories. For example, it is believed that 3d 
mirror symmetry arises when one rewrites three-dimensional supersymmet- 
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ric QED in terms of local operators which create Abrikosov-Nielsen-Olesen 
vortices |10|. This means that such operators are monopoles. However, it 



was never clear how to define monopole operators in SUSY QED. A proposal 
in this direction was made by M. J. Strassler and one of the authors |TI| , but 
it was only partially successful. 

In this paper we will address the second issue in a simple toy model: 
three-dimensional QED with Nf flavors of fermions. This theory is believed 
to flow to an interacting conformal fixed point for large enough Nf (this is 
discussed in more detail in the next section). The theory is not supersym- 
metric and is not expected to possess a simple dual. Nevertheless, we believe 
it is a useful exercise to define monopole operators in this simple model and 
learn to work with them. Besides, monopole operators are rather interesting 
beasts even in the non-supersymmetric case. First of all, these are the first 
examples of local operators in a three-dimensional CFT which are not poly- 
nomial in the fundamental fields. Thus our construction can be regarded 
as a generalization of the vertex operator construction from free 2d CFT 
to an interacting 3d CFT. Second, we show that because of fermionic zero 
modes monopole operators transform in a non-trivial representation of the 
flavor group whose size depends on the Chern-Simons coupling. Monopole 
operators in supersymmetric QED and their role in mirror symmetry will be 
discussed in a forthcoming publication. 



2 Review of three-dimensional QED 

The action of three-dimensional QED in the Euclidean space is given by 

Lqed = J d 3 x (J^F^ + $ (a • iD A ) ifA . 

Here A is the U(l) gauge field, F = dA is the field-strength 2-form, Da is 
the corresponding covariant derivative, and ip^ is a complex two-component 
spinor. The index j runs from 1 to Nf. 

In three dimensions one can add to the action a Chern-Simons term 

Such a term breaks parity invariance of the theory. We will assume that the 
gauge group is compact (i.e. U(l) rather than R). Naively, this requires the 
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Chern-Simons coupling k to be an integer, to avoid global anomalies. The 
real story is slightly more complicated. When Nf is odd, the fermionic path 
integral is anomalous. The anomaly is the same as the anomaly due to a 
Chern-Simons term with k = 1/2. Thus cancellation of global anomalies 
requires 

k ^- G Z. 

In particular, for odd Nf the Chern-Simons coupling must be non-zero, and 
parity is broken. This is known as parity anomaly ]12 . 



The gauge coupling e has dimension m 1//2 . Thus the theory is super- 
renormalizable and free in the ultraviolet. In fact, its UV behavior is so 
good that no renormalization of the Lagrangian is required. Contrariwise, 
3d QED is strongly coupled in the infrared (perturbative expansion is really 
an expansion in powers of e 2 /p, where p is the momentum scale). It is natural 
to assume that the low-energy limit of this theory is a non-trivial CFT, but 
this has not been conclusively demonstrated. However, the statement holds 
to all orders in 1/Nf expansion |T3], [TJ], |T3], |TB[. In fact, in the limit Nf — > oo 
the infrared theory becomes weakly coupled, and conformal dimensions of 
all fields can be computed order by order in 1/Nf. For example, the IR 
dimension of t/> 3 ' is canonical (i.e. the same as the UV dimension) up to 
corrections of order 1/Nf. 

More interestingly, the IR dimension of F^ v is 2 to all orders in 1/Nf. To 
understand why this is the case, consider a current 

= e ^PF vp . 

It is identically conserved by virtue of the Bianchi identity. A priori, this 
current could either be a primary field, or a descendant of the primary field. 
In the UV, the latter possibility is realized, since we can write 

J" = d»<r, (1) 

where a is a free scalar field. The scalar a is usually referred to as the dual 
photon. It has dimension 1/2 (as befits a free scalar in three dimensions), 
while and F^ v have dimension 3/2. On the other hand, in the IR an 
equation like Eq. (P cannot hold. Indeed, Eq. ([]]) implies that F pu obeys 
the free Maxwell equation, which clashes with the assumption that there are 
massless charged particles in the infrared. (We assume here that the fermions 
do not get a mass due to some non-perturbative effect, see a discussion 
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below.) This strongly suggests that in the IR limit J M is a primary field. 
It is well known that in a unitary 3d CFT a conserved primary current has 
dimension 2. Hence the IR dimension of J M and is 2. This conclusion can 
also be reached by directly studying the perturbative expansion in powers of 

Note that the difference between the UV and IR dimensions of F is of 
order 1, and therefore the IR fixed point is far from the UV fixed point, even 
in the limit Nf — > oo. In this respect, the situation is very different from the 



Banks-Zaks-type theories in four- dimensions [[17]], where the IR dimensions 
of all operators are very close to their UV dimensions. 

The physics of 3d QED at finite Nf remains controversial. The conven- 
tional approach is to study the system of Schwinger-Dyson equations trun- 
cated in some way and look for symmetry-breaking solutions. For simplicity, 
let us focus on the case of zero Chern-Simons coupling and even Nf. It 
has been claimed that at finite Nf flavor symmetry and parity are spon- 
taneously broken by a dynamical mass for the fermions, and the infrared 
limit is a free theory of photons ||18|| . The majority of such studies indicate 
that this happens for Nf smaller than a certain critical value of order 6 or 
7 (see e.g. J19], |20], |2T|). There are also claims that dynamical mass gen- 
eration takes place for all Nf but is exponentially small for large Nf and 
therefore invisible in l/Nf expansion [T^, |22|, It must be stressed that 
the results of such studies depend on the way one truncates an infinite system 
of Schwinger-Dyson equations, a procedure which cannot be fully justified. 
Lattice simulations of 3d QED have been inconclusive so far. 

In this paper we will be interested in the large Nf limit, and therefore the 
behavior at finite Nf will be unimportant. Note also that in the N = 2 and 
N = 4 supersymmetric cases the situation is better, in the sense that one 
can argue the existence of a non-trivial CFT at the origin of the quantum 
moduli space for all Nf. 
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3 Defining monopole operators 



3.1 A preliminary definition 

As mentioned above, three-dimensional QED possesses an interesting con- 
served current, the dual of the field strength: 

m — _L^p p 
"4tt vp - 

Its conservation is equivalent to the Bianchi identity dF = 0. The corre- 
sponding charge is called the vortex charge, because in the Higgs phase it is 
carried by the Abrikosov-Nielsen-Olesen (ANO) vortices. The vortex charge 
is integral if the gauge field A is a well-defined connection on a U(l) principal 
bundle. Loosely speaking, we would like to construct a vortex-creating oper- 
ator. But in an interacting conformal field theory, it does not make sense to 
say that an operator is creating a particle. A vortex-creating operator will 
be defined as an operator with a unit vortex charge. This means that the 
OPE of such an operator with J M has the form 

1 

J^(x)O(0) ~ — i — prO(0) + less singular terms. 

4:71 \X\ 6 

Such operators can be organized in the representations of the conformal 
group. In a unitary theory local operators must transform according to 
lowest-weight representations, i.e. those representations in which the di- 
mension of operators is bounded from below. The operator with the lowest 
dimension is called a conformal primary. It is standard to label a repre- 
sentation by the spin and dimension of its primary. Our problem can be 
formulated as follows: determine the spin, dimension, and other quantum 
numbers of primaries with a given vortex charge. 

In the path integral language, an insertion of an operator with vortex 
charge n at a point p is equivalent to integrating over gauge fields which 
have a singularity at x = p such that the magnetic flux through a 2-sphere 
surrounding x = p is n. To be consistent, one must regard charged matter 
fields as sections of a non-trivial line bundle on the punctured 1R 3 . Thus an 
insertion of a vortex-creating operator causes a change in the topology of 
fields near the insertion point. In what follows we will use the terms "vortex- 
creating operator" and "monopole operator" interchangeably. 

This way of defining topological disorder operators is familiar from 2d 
CFT. For example, a twist operator for a free fermion in 2d is defined by the 
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condition that the fermion field changes sign as one goes around the insertion 



point ||24|| . Another example is afforded by the theory of a periodic free boson 
in 2d. This theory has winding states, and the corresponding operators create 
a logarithmic singularity for the boson field. Thus our monopole operators 
can be regarded as three-dimensional analogues of twist operators or winding- 
state operators. 

In the two-dimensional case one can loosely say that a winding-state op- 
erator creates a kink. The precise meaning of this statement is the following. 
Consider a perturbation of the free boson theory by a periodic potential, say, 
a sine-Gordon potential. The resulting massive theory has multiple vacua 
and topological excitations (kinks) interpolating between neighboring vacua. 
The operator which carries winding number one has non-zero matrix elements 
between the vacuum and the one-kink state. 

Similarly, one can loosely say that a monopole operator creates an ANO 
vortex. To make this statement precise, one has to go to the Higgs phase 
(for example, by adding charged scalars with an appropriate potential). In 
the Higgs phase, the magnetic flux emanating from the insertion point of the 
monopole operator is squeezed into a thin tube. This tube is the world-line 
of a vortex. 



3.2 A more precise definition 

The definition of monopole operators given above is not yet complete. In 
effect, we have defined an insertion of a monopole operator by requiring that 
the gauge field strength have a particular singularity at the insertion point. 
However, we did not specify the behavior of the matter fields near the inser- 
tion point. In fact, we expect that there are many operators which carry the 
same vortex charge, and they differ precisely by the behavior of fields at the 
insertion point. Another difficulty is that the IR theory is strongly coupled, 
and it seems hard to compute correlators involving monopole operators. 

The first difficulty can be circumvented using radial quantization. It is a 
general feature of CFT in any dimension that local operators are in one-to- 
one correspondence with states in the Hilbert space of the radially quantized 
theory. This follows from the fact that one can use a conformal transforma- 
tion to map an insertion point of an operator to infinity. In this way one 
trades a local operator for an incoming or out-going state. In the case of 
monopole operators, such a mapping takes an operator with vortex charge n 
to a state on S 2 x IS with a magnetic flux n through S 2 . Here K. is regarded 
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as the time direction. Classifying states of a CFT on S^xR with a given 
vortex charge is certainly a well-defined problem. Furthermore, the radially 
quantized picture is the most convenient one for computing correlators which 
involve two monopole operators with opposite vortex charges and an arbi- 
trary number of ordinary operators. By mapping the insertion of a monopole 
operator to an in-going state and the insertion of an anti-monopole operator 
to an out-going state, one reduces the problem to computing a particular 
matrix element of a product of several ordinary operators. A particularly 
important special case is the three-point function which involves a monopole 
operator, an anti-monopole operator, and a conserved current. Knowledge of 
such correlators allows one to read off the quantum numbers of a monopole 
operator. For example, in order to determine the dimension of an operator, 
one has to compute the expectation value of the stress-energy tensor in the 
corresponding state. This approach is familiar from 2d CFT, where it is used 
to compute the quantum numbers of twist operators (see e.g. Ref. p4fl ). 

Of course, if one desires to compute four-point functions of monopole op- 
erators, mapping two of the insertion points to infinity does not help very 
much. In the case of 2d CFT, one has to use tricks special to the theory in 
question in order to compute four-point functions of topological disorder op- 
erators. In this paper, we will be content with studying the OPE of monopole 
operators with conserved currents, and leave the study of four-point functions 
to future work. 

The second difficulty can be avoided by working in the large Nf limit. It is 
a general feature of this limit that the gauge field does not fluctuate, and can 
be treated classically |13|, [14], [16| . This can be seen as follows. The infrared 



limit in 3d QED is simply the limit e — > oo. This is literally true, because no 
renormalization of the Lagrangian is required. Thus one can simply drop the 
kinetic term for the gauge field. Integrating out the fermions then gives an 
effective action for the gauge field of order Nf. For example, when expanded 
around a trivial background, this action looks like 

Nf J (F^ u U~y 2 F^ v + higher-order terms) d 3 x. 

Thus the effective Planck constant is of order 1/Nf, and in the large Nf limit 
the size of gauge-field fluctuations is order 1/Nf. Moreover, if we absorb a 

1/2 

factor of Nj- into F, we see that self-interactions of F are suppressed in the 

1/2 

large Nf limit. In other words, NJ F is a Gaussian field in the large Nf 
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limit. It is this line of reasoning that allows one to show that the infrared 
CFT is weakly coupled in the large Nf limit. The argument also applies to 
CFT on S 2 x R with a flux. Thus we can regard the gauge field as a classical 
background. 

It is very plausible that the saddle point of the effective action for F on 
S 2 x R is rotationally symmetric. Therefore we can assume that the classical 
background is simply a constant magnetic flux on S 2 . 

The above discussion reduced our problem to studying free fermions on 
S 2 x M in the presence of a constant magnetic flux. This is almost a textbook 
problem, and everything of interest can be computed. For example, finding 
the dimension of a monopole operator is equivalent to computing the Casimir 
energy of free fermions on S 2 with a flux. It is a priori clear that this energy 
scales like Nf. There are corrections to this result, which can be computed by 
taking into account the fluctuations of the gauge field. However, such effects 
are suppressed by powers of l/Nf. 

The above discussion contains a gap as regards gauge invariance of monopole 
operators. Gauge-invariance of a local operator is equivalent to gauge-invariance 
of the corresponding state in the radially quantized picture. In other words, 
the state must satisfy the Gauss' law. Gauss' law in QED on S 2 x R comes 
from varying the action with respect to the "time-like" component of the 
gauge field A. In the limit e — > oo it simply reads 



is the electric charge density operator. In particular, the total electric charge 
of a gauge-invariant state must be zero. The latter is a standard fact about 
gauge theory on a compact space, valid irrespective of the value of e. The 
definition of the electric charge operator involves normal-ordering ambigu- 
ities, which will be dealt with below. Note also that the inclusion of the 
Chern-Simons term in the action modifies the Gauss' law constraint into 



In particular, the total electric charge of the matter modes must be equal to 
—k times the vortex charge. In this way (and only in this way) the Chern- 
Simons term will affect the physics at large Nf. 



p(x)\$) = 0, 



where 
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4 Properties of monopole operators 



4.1 Radial quantization in the presence of a flux 

As explained in the previous section, at large Nf all properties of monopole 
operators can be deduced from the properties of free fermions on S 2 x R in 
a constant background magnetic flux. In this subsection we summarize the 
physics of this system, with detailed derivations relegated to the Appendix. 

The spectrum of the Dirac Hamiltonian on S 2 x K. with n units of magnetic 
flux is given by 

E p = ±y/p 2 +p\n\, p = 0,1,2,.... 

The degeneracy of the p-th eigenvalue is 2j p + 1, where 

1 

Jp = T^N ~ 1 )+P- 

These 2j p + l states transform as an irreducible representation of the rotation 
group SU(2) rot . 

The presence of n states with zero energy is particularly important. The 
existence of at least n zero modes is dictated by the Atiyah-Singer index 
theorem applied to the Dirac operator on S 2 coupled to the magnetic field. 

In the case of a unit magnetic flux (|n| = 1), we have a single fermionic 
zero mode with zero spin. Thus a spinor is converted into a scalar due to 
the non-trivial topology of the magnetic monopole. This scalar-spinor trans- 
mutation is well known in other contexts; in particular it plays an important 
role in the conjectured S-duality of iV = 4 d = 4 super- Yang-Mills. For gen- 
eral n, the fermionic zero modes transform in an irreducible representation 
of SU(2) rot with spin j = (\n\ — l)/2. We will discuss in detail the case when 
n = ±1, and then comment on the higher-n case. 

Let us denote the fermionic annihilation operators by c pm , where i = 
\,...,Nf is the flavor index, p = 1,2,..., labels the energy eigenspaces as 
above, and m = —j p , —j p + 1, . . . , j p , labels the states within the p-th energy 
eigenspace. The fermion annihilation operators corresponding to p = will 
be denoted simply by c z Q . The Hilbert space of the theory is the tensor product 
of the Hilbert space of zero modes and the Hilbert space of all other modes. 
The latter is simply a fermionic Fock space with a unique vacuum \vac) + 
which satisfies 

c pm \vac) + = 0, j9>0,Vi,m. 
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This vacuum state is rotationally invariant. 

The Hilbert space of zero modes is also a Fock space of dimension 2 N f , 
with the vacuum vector which we denote \vac)o. It is spanned by the vectors 

\vac)o, c^\vac) , Cq> \vac) , c^c^ . . . c \vac} . 

All these states are degenerate in energy, and none is a preferred vacuum. 
Since the zero modes have spin zero, all the ground states are rotationally 
invariant. We conclude that the radially-quantized theory of free fermions 
has a 2 N f -fold degenerate ground state. 

However, we still need to impose the Gauss' law. The charge density 
operator receives contributions from both zero and non-zero modes. The 
part due to non-zero modes can be defined using the obvious normal-ordering 
prescription. If we put all non-zero modes in the vacuum state, then the 
charge density due to non-zero modes vanishes. It remains to analyze the 
contribution from zero modes. Naively, it seems that the Fock vacuum must 
be assigned zero electric charge. If this were the case, then the states obtained 
by acting on the vacuum with zero mode creation operators would have 
positive charge, and therefore would not be gauge-invariant. But because of 
normal-ordering ambiguities, the situation is more interesting. 

As stressed above, the Fock vacuum for the zero modes is not that special. 
The completely filled state appears to be an equally good candidate for a 
state with vanishing electric charge. The two just differ by a change in 
the normal ordering prescription. A statement which is independent of the 
normal-ordering prescription is that the electric charge of the filled state 
exceeds the charge of the vacuum by Nf. If one wants to be "democratic", 
one has to assign charge —\Nj to the vacuum and charge |iV/ to the filled 
state. A similar symmetric charge assignment has been advocated by Jackiw 
and Rebbi in their pioneering study of fermions bound to solitons, on the 
grounds on charge-conjugation symmetry [[25[ . 



The precise argument for the symmetric charge assignment goes as fol- 
lows. Charge conjugation maps a monopole to an anti-monopole and by 
itself does not tell us anything. But CP transformation maps a monopole 
to itself. If we want to quantize in a CP-invariant manner, we must assign 
opposite electric charges to states related by CP. Since CP takes annihilation 
operators into creation operators, the filled state and the vacuum are related 
by CP, and their electric charges must be opposite. 

The invocation of CP invariance assumes that the theory we started with 
is CP-invariant. This means that the symmetric charge assignment is valid 
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for a vanishing Chern-Simons coupling. But we know that turning on the 
Chern-Simons coupling k is equivalent to shifting the electric charge by k 
times the vortex charge. Therefore we conclude that in the presence of the 
Chern-Simons coupling the Fock vacuum has electric charge 



while the filled state has charge 



Note that because of the parity anomaly, the electric charge is all always 
integer-valued, whether Nf is even or odd. This a manifestation of the close 
relationship between the existence of parity anomaly and the induced vacuum 
charge |p6| . 

Now we can analyze the consequences of the Gauss' law constraint. If all 
non-zero modes are in their ground state, the constraint simply says that the 
total electric charge of the state must be zero. For k = this implies that a 
physical state is obtained by acting with Nf/2 zero modes on the vacuum. 
The number of such states is 



and they transform as an anti-symmetric tensor oiSU(Nf) with Nf/2 indices. 
Note that cancellation of global anomalies requires Nf to be even when k = 0, 
so this result makes sense. For k between —Nf/2 and Nf/2 the physical 
states are obtained by acting with Nf/2 — k zero modes on the vacuum. The 
corresponding states transform as an anti-symmetric tensor of SU(Nf) with 
Nf/2 — k indices. Again global anomaly cancellation ensures that Nf/2 — k is 
an integer. For \k\ > -y- there are no gauge- invariant states with unit vortex 
charge and all non-zero modes in their ground state. 

If one does not assume that positive-energy modes are in their ground 
state, then one can construct many other states which satisfy the Gauss' law 
and have unit vortex charge. However, such states will have higher energy 
than the ones discussed above. 

Now let us consider the more complicated case of n — 2. For simplicity 
we will set the Chern-Simons coupling to zero and take Nf to be even. In 
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the case n = 2 each fermion has two zero modes which transform as a spin-| 
representation of SU (2) rot . Reasoning based on CP-invariance tell us that the 
Fock vacuum has electric charge — Nf. Physical states must have zero electric 
charge and are obtained by acting with Nf zero modes (out of a total number 
of 2Nf) on the vacuum. But physical states must also be annihilated by the 
electric charge density operator. This is not automatic anymore, because the 
fermionic zero modes are not rotationally invariant. A short computation 
shows that the electric charge density operator for the zero modes po(x) has 
a piece which transforms as a singlet of SU (2) rat and a piece which transforms 
as a triplet of SU(2) rot . The former is simply the average of po(x) over the 
sphere and is proportional to the total electric charge. The spin-triplet piece 
of po(x) is proportional to the total spin, simply because this is the only 
spin-triplet one can make out of two spin-1/2 fermions. Thus the Gauss' law 
constraint is equivalent to the requirement that the total electric charge as 
well as the total spin be zero. 

For example, for Nf = 2, there are six states with zero total electric 
charge, which are obtained by acting on the Fock vacuum with two zero 
modes out of the available four. Three of these states transform as a vector 
of SU (2) rot and as singlets of the flavor group SU (2) f and do not satisfy the 
Gauss' law constraint. The remaining three transform as singlets of SU(2) rot 
and as a triplet of SU(2) f. These three states are gauge-invariant. Note that 
in this case gauge-invariant states transform as an irreducible representation 
of the flavor group. For Nf > 2 this is no longer true, as one can easily check. 

4.2 Quantum numbers of the monopole operators 

In this section we determine the quantum numbers of the simplest monopole 
operators, the ones with the lowest conformal dimension for a given vortex 
charge. On general grounds, such an operator lives in a lowest-weight rep- 
resentation of the conformal group, and its conformal dimension is defined 
as the conformal dimension of the lowest-weight vector, or, if we pass to the 
radially quantized picture, as the energy of the corresponding state. 

Let us begin with the case n — 1. As explained above, gauge- invariant 
states with lowest energy are obtained by putting all non-zero modes in their 
ground states and acting by Nf/2 — k zero mode creation operators on the 
vacuum. Obviously such states transform as an anti-symmetric representa- 
tion of SU(Nf) with Nf/2 — k indices. It is interesting to note that the 
usual gauge-invariant operators which are polynomials in the fundamental 
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fields transform trivially under the center of SU(Nf). Indeed, free fermions 
have flavor symmetry group U(Nf), and since we are gauging its U(l) sub- 
group, the flavor symmetry of QED appears to be U(Nf)/U(l) = PU(Nf) = 
SU(Nf)/%Nt- But monopole operators transform non-trivially under Zjvy 
(except for k = ±iV//2). A very similar effect occurs in N = 2 d = 4 super- 
QCD, where all perturbative states transform as tensor representations of 
the flavor group SO(2Nf), while magnetically charge states transform as 
spinors 
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Other quantum numbers of interest are spin and conformal dimension. 
Since the Fock vacuum and the zero modes are rotationally invariant, the 
spin of our monopole operator is zero. The dimension is proportional to the 
energy of the state. As usual, the definition of the energy is plagued by 
ordering ambiguities. However, we have a simple cure: we can normalize the 
energy by requiring that the unit operator have zero dimension. This means 
that the energy of the ground state on S 2 with zero magnetic flux is defined 
to be zero. The energy of any other state can be defined by introducing a 
UV regulator, subtracting the regularized energy of the state corresponding 
to the unit operator, and then removing the regulator. This procedure gives 
a finite answer, which is not sensitive to the precise choice of the regulator, 
provided the regulator preserves the symmetries of the problem. 

In order to make precise the relation between the Casimir energy and the 
dimension, recall that the OPE of a spin-zero primary field and the stress- 
energy tensor reads: 



T. u (x)0(y) ~ — ( J 9 d 1 \ Q , n 



where h is the conformal dimension. If the stress-energy tensor of free 
fermions is defined by 

% — % — — 

= -~i/> (j^V v + 7„D M ) ^ + - (£^7^ + T>v$l v ) $ ~ 9»v£, 

then h^p — h$ — 1, the standard normalization. This implies that in the 
radially-quantized picture the expectation value of the stress-energy tensor 
in the state \0) is given by 



{T. v dx» <g> dx u ) = ^-(dT 2 - -{d6 2 + sind 2 
An \ 2 
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Thus h is simply the energy of \0) with respect to the Killing vector J^. In 
our case, this means that the conformal dimension of the monopole operator 
is the Casimir energy of Nf free fermions on S 2 with a magnetic flux. This 
Casimir energy for any n is computed in the Appendix. For n = 1 the result 
is 

h 1 = N r 0.265 .... 

By charge-conjugation symmetry, monopole operator with n — — 1 has the 
same conformal dimension and spin and transforms in the conjugate repre- 
sentation of the flavor group SU(Nf). 

It is easy to extend the discussion to n — ±2. As explained in the 
previous section, Gauss' law constraint is equivalent to the requirement of 
zero spin and zero electric charge. The states with zero electric charge are 
obtained by acting with Nf zero modes (out of total number of 2Nf zero 
modes) on the Fock vacuum. These states transform as an anti-symmetric 
tensor of SU(2Nf) with Nf indices. Gauge- invariant states are obtained 
by decomposing this representation with respect to the SU(2) rot x SU(Nf) 
subgroup and separating out S , f/(2) ror singlets. In general, gauge-invariant 
states transform as a reducible representation of SU(Nf). One can easily 
show that the dimension of this reducible representation is 

( \N} + N f -1\ 
V 5*/ J 

The conformal dimension of the corresponding monopole operators is the 
Casimir energy of free fermions in a background magnetic field. Numerically, 
it is given by 

h 2 = N r 0.673 .... 

It is interesting to note that 2h\ < h 2 (at least for large Nf). Therefore 
the OPE of two monopole operators with n = 1 and the lowest conformal 
dimension contains only terms with positive powers of \xi — x 2 \. 

5 Discussion 

In this paper we have constructed local operators in an interacting 3d CFT 
which carry vortex charge and therefore create Abrikosov-Nielsen-Olesen vor- 
tices in the Higgs phase. We have shown that for large Nf such operators 
have conformal dimensions of order Nf. For the case of unit vortex charge, 
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we showed that the operator with the lowest possible dimension has zero 
spin and transforms in a non-trivial representation of the flavor group. An 
important tool in this analysis is 1/Nf expansion. 

The idea that vortex-creation operators can be studied in the large Nf 
limit has been previously proposed in Ref . |Tl| . The approach taken there was 
to integrate out the matter fields, and then perform a duality transformation 
on the effective action for the gauge field. Then the vortex-creation operator 
is defined as the exponential of the dual photon. One drawback of this 
approach is that it is easy to miss fermionic zero modes, and consequently 
to misidentify the quantum numbers of the vortex-creating operator. It is 
preferable to keep the matter fields, and to identify a vortex-creating operator 
by the property that its insertion causes a change in the topology of the gauge 
field. As we have seen above, this definition can be made concrete by using 
radial quantization and large Nf expansion. 

Our main motivation for studying vortex-creation operators was the hope 
that this would enable us to give a constructive proof of 3d mirror symmetry. 
It is straightforward to apply the methods of this paper to 3d gauge theories 
with N = 2 and N = 4 supersymmetry. The results and their implications 
for mirror symmetry will be reported in a forthcoming publication. 

It is natural to wonder if our approach to the construction of topological 
disorder operators has an analogue in four dimensions. In three dimensions, 
we defined the vortex charge of a local operator as the first Chern class of 
the gauge bundle evaluated on an S 2 surrounding the insertion point. In 
four dimensions, we have S 3 instead of S 2 , and since characteristic classes of 
vector bundles are even-dimensional, it appears impossible to define a similar 
topological charge for local operators. On the other hand, a B-field on an 
S 3 can have non-trivial topology, since its field-strength is a 3-form. Thus, if 
there were an interacting 4d CFT involving a B-field, one could define local 
operators which create topological disorder. In order for this to work, the 
field-strength 3-form must have dimension 3, so that its dual is a conserved 
primary current. Note that in the theory of a free B-field, the field-strength 
has dimension 2. In this case the dual current, although conserved, is not a 
primary, but a gradient of a free scalar. Thus in order to define a conformally- 
invariant topological charge, the 4d CFT must be interacting. Unfortunately, 
no such theory is known at present. Perhaps there exists a duality-symmetric 
reformulation of N = 4 d = 4 super- Yang-Mills which involves B-fields, and 
in which both W-bosons and dual W-bosons are described by topological 
disorder operators. 
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After the first version of this paper was posted on the arXive, we learned 
that topological disorder operators in 3d have been previously considered by 
G. Murthy and S. Sachdev p8"] . The model considered there was the CP^ 
model in three dimensions. This theory is not renormalizable and requires an 
ultraviolet cut-off. It has a critical point separating the ordered phase, where 
the sigma-model field has an expectation value, and the disordered phase, 
where the correlators decay exponentially. Unlike in 3d QED, there are dy- 
namical topological defects in the CP^ model (so-called hedgehogs). But at 
the critical point and in the limit of large N their density vanishes, and the 
situation becomes very similar to that in 3d QED. In particular, the critical 
exponents computed in Ref. ]28| can be interpreted as scaling dimensions of 
hedgehog operators. It is interesting to note that the approach to computing 
these scaling dimensions taken by Murthy and Sachdev is rather different 
from ours. Instead of evaluating the expectation value of the stress-energy 
tensor in the presence of topological disorder operators, they in essence com- 
pute the 2-point function of these operators. This is somewhat obscured by 
the fact that Murthy and Sachdev map both operator insertions to infinity. 
As a result, the distance dependence of the 2-point function is traded for an 
anomalous dependence of the partition function on the ratio of the ultravi- 
olet cut-off A and the infrared cut-off A. This method could be used in 3d 
QED as well. In fact, M. J. Strassler and one of the authors of the present 
paper (A.K.) have contemplated such a route to computing scaling dimen- 
sions of monopole operators in 3d QED and SQED, but were discouraged by 
apparent technical difficulties. It would be interesting to rederive the results 



of the present paper using the approach of Ref. [28 



Appendix 
Monopole harmonics 

To solve for the energy spectrum of free fermions on S 2 with a magnetic flux, 
we will use the fact that this system is related by a conformal transformation 
to the Dirac equation on M 3 in the monopole background. This allows us to 
use the machinery of "monopole harmonics" developed by Wu and Yang [23 
The three-dimensional Dirac operator on flat M. 3 is given by 

iD = —a • ff, 
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where cr. 



a y , and a z are the Pauli matrices, and 7? 



p + A, with p being 



the momentum operator. Following Ref. p9| , let us define the generalized 
orbital angular momentum operator as 



L = r x 7r 



with q = —eg = n/2. It is straightforward to check that L defined this way 
satisfies the angular momentum algebra: 



[Lj , Xk] 

[Lj,lTk] 

[Lj,L k ] 



We define the total angular momentum as 



J 



We can take r, L 2 , J 2 , and J z as a complete set of observables (it is easy 
to check that they commute and are all self-adjoint with respect to the usual 



J AD 



0, but 



L , ID 



^ 0. How- 



inner product). It can be checked that 

ever, this is good enough, because as we will see later, to find the eigenvalues 
of %D we only need to diagonalize an operator in a two-dimensional space. 
The monopole harmonics Y q ^ m (Q, 4>) constructed in Ref. p{| satisfy 



1 = \q\, \q\ + i, \q\ + 



2,..., 



I J zYq,l,m 

m = 



mY Q i 



q,L,mi 
-/,...,/. 



The simultaneous eigenfunctions of |L , J , J z j will be denoted by 4>ij mj and 
are given by 




l—m XT 
2i + i I q,l,m 

'l+m+ly 

21+1 1^,^+1 



l+m+l y ■ 

tor j = I H — (m,- = m H — J 
_jhv , , , / 2 2 



1 1 

for j = I - - (I ^ 0, rrij = m + -). 
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They satisfy 

L 2 4>lj mj = 1(1 + l)(f>ljmj, 
J 2 (Pljm, = j(j + l)</>lj mj , 

We can summarize the possible value of l,j,rrij as follows: 

1 . 1 3 5 

• ' = M-2' |g| + 2' |g| + 2' M + 

(for q — 0, j = |g| is not allowed); 

Z 

1 1 1 

• if j — \q\ — -, then / = j + - = \q\, otherwise I — j ± -; 

Any wave-function can be expanded as 

= R ijm,(r)(f)i jmj (e,(f)). 

l,j,mj 

Note that while <f)ij mj is a two-component spinor, Rij mj is just a scalar. 
Now we may write iD in terms of the angular momenta. Define ay as 

a • r 

r 

One can show that 

d 1 _ =f a r 

aJiD) = i— i-a • L — %q — , 

or r r 

where we made use of the fact that 

(a ■ G) (a ■ K^j = G ■ K + ia ■ (d x it} 

for any G and K that commute with a. Now using the fact that oy 2 = 1, we 
have 

•n fr>\ ■ 9 - Gr ^ t ■ 1 
iD = o r o r [%D) = ia r — % — cr ■ L — iq— 

or r r 

= ia r - i — (J -L --\-iq-. 

Or r 4 r 
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Thus the Dirac Lagrangian on R 3 in the presence of a monopole can be 
written as 

C Ri = % -i)(j r - {J 2 ~ L 2 - |) - qa^j ip. 

Setting r = e T and performing a Weyl rescaling 

we obtain the Lagrangian on S 2 x R: 

£ S 2 X r = #<r r ( - (J 5 - L 2 + - gov J ip. 



Note that the norm 



L 



r 2 dflipa r ip 



S2 

on R 3 is transformed to the norm 



/ IpCTi 
JS2 



on S 2 x R. 

Taking into account the above results, the Euclidean equation of motion 
for ip is 



dRij mj (r) 
dr 



- + + + -)R ljmj (r) 



- ^ qRi'j'mri^iljrrijlarll'j'm'j) = 0, 



I'j'm'- 



where we have used {ljmj\a r \l'j'm'j) to denote J dQ<p\- m .a r <^i> j' m '.. 



Now the identity 



</, <J r 



tells us that 



(Ijrrijlarll'j'm'j) = 6 jj >6 m . m '.(ljm j \<T r \l'jm j ), 
and thus the eigenvalue equation becomes, for any given j,m,j, 

dRljmj (t) 

~dr 



+ -1(1 + 1) + -)R ljmj (r) 



- ^ qRyjm^^jm^a^l' jrrij) = 0. 
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Let us suppress the j, rrij indices, and denote R(i=j-i)j mj by R a , \l — j — 
\jrrij) by \a), R {l=j+ i )jm . by R h , \l = j + ^jnij) by |6>, (a\a r \a) by a aa , 
(a\a r \b) by a a b, (b\a r \a) by a ba , and (b\a r \b) by a b b- Then for any given j,m,j, 
we have two coupled first-order differential equations: 

<W " (T) - fj + l)R^r) + q (a aa R a (T) + a ab R b (r)), 



dr V 2 
dR b (r) ( . 1 



- (j + §J i? 6 (r) + g(a bbJ R b (r) + a 6a i? a (r)). 

A straightforward calculation of the matrix elements a aa , <r a b, and cr^ 
gives 

-Q 1 L ( Q 

J + 2 ^ + 2 V VJ + 2 

and of course a ba = cr* b = a ab . 

Energy spectrum and Casimir energy for fermions on 

S 2 

The energy spectrum can be read off from the behavior of the solutions as a 
function of r: a solution with energy E behaves as e~ Er . The results are as 
follows. 

Case (i): q = 0. 

The two equations decouple, and we find 

R a ( T ) = C a e^\ R\r) = C b e^ j+ ^ T , 

where C a and C b are integration constants, and j = |, |, | . . . There are no 

zero-energy solutions. 

Case (ii): q ^ 0, j = \q\ — |. 

In this case, there is no such thing as R a (because I cannot be j — | = 
\q\ — 1). So the first equation is absent, and the second equation gives: 

R\r) = C, 

with an arbitrary constant C . This solution has zero energy and degeneracy 
2j + l = 2|g|. 
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Case (Hi): q ^ and j = \q\ — | + p > \q\ — \, p — 1, 2, 

In this case, the two equations are coupled but easy to solve by eliminating 
one of the two unknowns. The result is 



R a ( T ) = qC ie T V ( j+ ^Y-i 2 + qC 2 e- T V (^Y-i\ 

' 1 



R"t) 



] + 



c 



+ 



— q + [ j + 



C 2 e 



2/ V 2 , 

where Ci and C2 are integration constants. The corresponding energies are 

- q 2 = ±\f2\q\p + p 2 , 



E r . 



±\ 



] + 



with degeneracies 2j + 1 = 2\q\ + 2p. Note that the spectrum is symmetric 
under q — > — q. 

The regularized Casimir energy is given by 



E reg ((3) = -J2( 2 P+ \n\)^p 2 + p\n\e'^ p2+pH . 

p=0 

We renormalize it by requiring that the Casimir energy of the vacuum with 
n = be zero. That is, we subtract from the above sum a similar sum with 
n — 0, and then take the limit (5 — > 0. Using the Abel-Plana summation 
formula 

F(it) - F(-it) 



J2f( P ) = -F(0)+ dxF{x) + i 
we obtain a finite answer for the Casimir energy: 

Ecasimir = + M(M - 2) + 

^jf* ( !i + T + 1 )v 



- 1 



Inl \ 2 n 2 



= 27Tt 



- 1 



Here one needs to take the branch of the square root which is positive on the 
positive real axis. The integral cannot be expressed in terms of elementary 
functions, but can be easily evaluated numerically for any n. 
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